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c^ \ We derive quantum Boltzmann equations from Schwinger-Dyson equations in gradient ex- 

pansion for a weakly coupled scalar field theory with a spatially varying mass. We find that 
at higher order in gradients a full description of the system requires specifying not only an 
on shell distribution function but also a finite number of its derivatives, or equivalently its 
higher moments. These derivatives describe quantum coherence arising as a consequence of 
localization in position space. We then show that in the limit of frequent scatterings coherent 
quantum effects are suppressed, and the transport equations reduce to the single Boltzmann 
equation for particle density, in which particles flow along modified semiclassical trajectories 
in phase space. 



1 Introduction 

In this letter we present a controlled derivation of dynamical transport equations for a simple 
complex scalar theory 

£ = {d^<i))\d^<i)) - m\x, t)</>V + Ant , (1) 

where the mass represents coupling to a classical background field which varies in space and 
time, and £int denotes interactions. For definiteness we consider here the simple quartic 
interaction Ci^x. = ^-^(0*0)^/4. In particular we study the effect of the breakdown of trans- 
lational invariance with a treatment of the varying background mass at nontrivial order in 
gradients. 

The motivation for this work comes from electroweak baryogenesis at a first order phase 
transition [Q], where one needs to model the departure from thermal equilibrium induced 
at the phase boundary of a growing bubble of the broken phase. No systematic treatment 
of such plasma dynamics, adequate to incorporate the crucial CP violating effects, is yet 
available. This paper is the second in a series in which we attempt to provide a systematic 
formalism for the treatment of this problem. Including higher order gradients is required 
since, in realistic cases, the relevant CP violating effects occur typically at higher order in 
gradients. The formalism we are developing is general and may be applied to other problems 
which involve analogous physics. 

Our treatment of the problem is based on the out-of-equilibrium closed time contour 
(CTC) formalism. In our derivation we assume a weak coupling limit and the semi classical 
approximation. The semiclassical condition, fcL 3> 1, states that the de Broglie wave length 
in the relevant direction must be large in comparison to the corresponding scale of variation 
of the external field L. Since electroweak bubble walls are believed to be thick: L ~ (10 — 
20)T~^, and relevant particle species are weakly coupled, these constraints are satisfied for 
most of particles in the electroweak plasma. Our treatment does not apply to reflecting 
particles, for which A;L ~ 1. 

In a recent paper [Q we have analyzed the propagator of the scalar field theory in Eq. (0) 
to non-trivial order in gradient expansion. The main implication of this analysis is that the 
quasiparticle picture, in which the plasma is treated as a collection of one-particle excita- 



tions with a given dispersion relation, breaks down. The reason is that in the absence of 
translational invariance, states locahzed in coordinate space mix coherently in momentum 
space. As a consequence, a self-consistent description of transport of plasma excitations 
requires additional equations which encode information about quantum coherence. In this 
letter we derive such a set of quantum transport equations. Further we show that in the 
limit of frequent scatterings when TqcL ^ 1, where Tqc is the relevant scattering rate, the 
terms describing coherent quantum effects are suppressed and can be neglected. This can 
be understood simply as the result of the scattering projecting onto the local semiclassical 
(quasiparticle) states. In this limit we recover a Boltzmann equation which describes colli- 
sions and canonical flow on non-trivially modified semiclassical trajectories in phase space. 
In this letter we restrict our analysis to a simple scalar theory, but we expect that the struc- 
ture of the transport equations derived is generic. In the discussion section we briefly outline 
the relevance of our findings for computations of baryon production at the electroweak scale. 

2 Kadanoff-Baym Equations 

The basic quantity in our derivation is the 2-point Green function in the out-of-equilibrium 
field theory. It is most conveniently defined in the Schwinger-Keldysh formalism |^ on a 
closed time contour (CTC) 

Gcix,y) = -t{Tc[^ix)^\y)]), (2) 

where Tq defines time ordering along the contour C which starts at some to 5 often taken to 
be at —00, goes to +00, and then back to to- The two point function Gc{x,y) obeys the 
contour Schwinger-Dyson equation (see Fig. 1): 

Gc{x,y) = G'c{x,y) + J^dx' J^dx" G^c{^,x')Ec{x',x")Gc{x",y) , (3) 

where Ec is the self-energy and Gq is the free particle (tree level) propagator. In order to 
solve Gc{x, y) from (^ some external information about the self energy function S^ must be 
provided. In general this means coupling infinitely many new equations to (^, leading to the 
quantum generalization of the BBGKY hierarchy. In the weak coupling limit it is natural to 
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Figure 1: Diagrammatic representation of the Schwinger- Dyson equation for the scalar theory 
with C,^, = -\{<P*<Pf/A. 

truncate the hierarchy by substituting all higher than 2-point functions by the perturbative 
value for the corresponding interaction vertex. For the scalar theory in Eq. (|lD the simplest 
such truncation consists of approximating the four point function by the simplest four point 
vertex linear in the quartic coupling A. As a consequence the two-point self-energies are 
evaluated at two loops (c/. Fig 1) and the resulting dynamical equations are then truncated 
at the order A^ assuming the weak coupling limit. 

The complex time ordering formulation in (0-^ can be conveniently expressed in terms 
of the usual time ordered Green functions along the real axis. If we define 

G>{x,y) = -z(0(x)0(y)t) (4) 

G<{x,y) = -^(0(y)V(a;)) (5) 

G^Xx,y) = eixo-yo)iG>ix,y)-G<ix,y)) (6) 

G''{x,y) = -e{yo~Xo){G>{x,y)-G<ix,y)), (7) 

where G"^'" are the retarded and advanced functions, we can write (§) as an equivalent set of 
Kadanoff'-Baym eqnaXions [^: 

{G,'-J:^'^)^G^'%x,y) = S{x-y) (8) 

(Go-i-SO®G<'>(a;,2/) = S<'> ® G"(a;,y), (9) 

where Gq^ = d"^ — m? is the inverse of the free propagator and ® represents the convolution 
integral: A ® B{x,y) = J d'^zA{x, z)B{z,y). Functions (§-0) are not all independent. In 
fact they can be reduced to just two independent real functions by using (i) the hermiticity 
constraint G'''{x,y)* = G"'{y,x) and G^{x,y)* = -G<{y,x), (ii) the relation G'' - G"^ = 
G^ — G^ and (iii) the spectral representation which relates ImG^''^ to ReG"'''^. 

Transforming to the Wigner representation, G{k;X) = J d'^re^'^'^G{X + r/2,X — r/2), 
and making use of the relation (S ® G){k; X) = exp —iO{Tj{k; X)}{G{k; X)}, the real part 



of (H) becomes an equation for the propagator 

cos 0{fi2 ± iujT}{G'''''} = 1 (10) 

and the real part of the equation (^ yields the quantum Boltzmann equation for the dy- 
namical variable G"^: 

- sin 0{n^}{tG<}=- cos 0({S>}{G'<} - {S<}{G'>}) - sin 0{zS<}{G'^}, (11) 
where O is the Poisson bracket operator 

0{f}{9} = I [dxf ■ d,g - d^f ■ dxg] (12) 

and we defined the shorthand notation 

n"^ = k"^ - m^ - j:n. (13) 

The retarded and advanced operators were decomposed as 

G'■''^ = GRT^A, S'-'° = S^T*^r, (14) 

where A is the spectral function. Finally, assuming spectral decomposition of the Wigner 
functions and making use of ^ = i{G^ — G^)/2, we can write 

iG< = 2An iG> = 2A{n + 1). (15) 



Eqs. ([TO|-[TID are the full dynamical equations in the Wigner representation. These are 
suitable for the description of systems in slowly varying backgrounds, with truncation at a 
given order in gradients leading typically to a more accurate modeling of the dynamics. 

2.1 Propagator equation 

To the lowest order in gradients Eq. (p!0D defines the familiar propagators in a spatially 
constant background. At this order the spectral function A is singular in the limit F — i> 0, and 
defines the well known projection to the local quasiparticle on-shell. Given the decomposition 
(p!5|), the QBE ([TT|) also becomes singular, allowing a reduction of the QBE by integration 



over momentum (or frequency) to the well known semiclassical Boltzmann equation involving 
on-shell excitations only ||^. When higher order gradient corrections are included, the on- 
shell projection becomes more involved because of the coherent quantum effects described 
by the gradient terms [^. Before performing the on-shell projection for the QBE, we shall 
here illustrate the technique using a simple test function. 

Solving Eq. ([To| ) iteratively to the lowest nontrivial order in gradients around the lowest 
order pole gives 

^ z 2 z^ 2 z^ ' ^ ' 

where z = k^ — k^ with /cq = cj^ — ku — m?{x) and we have assumed S/j = for simplicity. 
The spectral integral over the momentum variable k of some function T can be converted 
to a contour integral encircling the multiple pole of G at 2 = 0: 
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Ip[T] = - dk^ApT 

TT JO 



Res 



p+i 



GpT/^ki-z =Y.c,r^\ko), (17) 



z=0 



i=0 



where T^^^ko) = {d\JT){kQ)^ and the coefficients Cj may contain gradients up to pth order 
with respect to x. To the leading order in gradients {p = 0) the spectral function indeed 
becomes singular: A -^ {'K/2ko)[6{kx — ko)+6{kx + ko)], projecting sharply on-shell k^ = ifco- 
It should be noted that the singularity of the propagator ( [I6| ) remains at the quasiparticle 
shell k1 = kQ. The effect of gradient corrections in Gp is to project out derivatives of T w.r.t. 
kx up to order p + 1. That this produces an effective shift of the pole emerges when the 
contribution from the first order derivative is included in the definition of the projected 
function. To second order in gradients 

UT] = ^^ + c,T^'\ko) + c,T^'\kol (18) 

where 

k^ =ko-\ --H ^ ,/ (19) 

" 8 ki 32 kl ^ ' 

is the space dependent semiclassical momentum which coincides with the standard WKB 
dispersion relation. The higher order derivative corrections in Eq. (|IBp, however, are of the 
same order in dx as the shift. As we will now see in detail this leads to the breakdown of the 
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quasi-particle approximation when the projection is performed on the QBE. It is then not 
sufficient to describe the system with a single distribution function obtained by projecting 
on-shell, but some off-shell information, represented by a finite number of derivatives of n, 
describing coherent quantum effects, is necessary. 

3 Quantum Boltzmann Equation 

We will now consider the QBE (0) to the lowest nontrivial order in gradients. We truncate 
the collision term at leading order in gradients and focus on the structure of the flow term in 
the presence of a varying background. This approximation amounts to neglecting the deriva- 
tives of the self-energies, while retaining those of the background. (A complete consideration 
that includes all second order terms in the collision term will be given elsewhere [^]). We 
then have 

-(o- -oA {n^}{iG<} = -[S>G< - S<G>]. (20) 

The on-shell projection of the quantum Boltzmann equation is usually done by integrating 
over frequencies, resulting in an equation on the phase space {k; t,x} 0. An alternative yet 
equivalent approach of integrating over momenta is more convenient here, because in the 
present spatially varying problem the density of states is most conveniently labeled by the 
conserved energy. 

Let us continue with the special case T,r = 0. Inserting the decomposition (|T3p into ( PU[ ) 
and writing the O-terms explicitly gives 

{^d, + k-d,+ ^{d,kl)d,^ - ^{dlkl)dl^^ An = -^A (zS>n - ^S<(n + 1)) . (21) 

Integrating (^) over k^ gives a dynamical equation coupling n{ko) and its ffist three deriva- 
tives. The information contained in this zeroth moment clearly does not suffice to define 
a closed solution to the problem. To provide closure we can perform integrals of Eq. ( pTD 
weighted by some higher powers of kx, in a manner analogous to the standard derivation of 
fluid equations by taking moments of the classical Boltzmann equation. There is however 
one crucial difference. While for the latter case there is no control parameter and hence no 
natural closure exists, in the former case the equations close at a finite number of indepen- 
dent moments. This is the case because, as shown in equation (0) above, in an integral over 
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any smooth test function weighted by Gp only the first p + 2 terms are nonzero. In particular 
at second order in gradient expansion the closure is obtained by the first four moments. 

The QBE becomes very complicated when written in terms of n^'^\kQ). It is more con- 
venient to express the equations in terms of the weighted projections of the generalized 
distribution function n 

f, = e{K)f^ + e{-h)fr, (22) 

where 

2 r°° 2 r^ 

/+ = - / dkXAn, fr = - dKi-KYAn. (23) 

These functions are a straightforward generalization of the distribution function for the 
spatially constant case. In the absence of gradient corrections f^ -^ kQ^^n(±ko), so that 
the standard distribution function is then simply OikQjniko) + 6'(— A;o)^(— ^o)- Now using 



moments (p3|) it is particularly easy to prove the closure. Indeed, at second order in gradients 



A ^ A2-, and we have the identity 

/ dk^{k^-ka)%A2n = Q, (24) 

Jo 

and the analogous identity holds for A;^: < 0, so that one immediately obtains 

/,+4 - 4A;o/i+3 + Qklfi+2 - ^k'ofi+i + ^tfi = 0- (25) 

The constraint (^) allows any moment // to be written in terms of the four lowest ones 
/o) /i; /2, /s- A similar binomial constraint holds at p-th order in gradients, providing closure 
with p + 2 moments. 

It is curious to observe that the functions fi can be interpreted as projections onto 
different momentum hypersurfaces ki of n. Performing the fc-integral in Eq. (pSf ) implies 
that 

/± = k[-^ Ki ni ±h) + 0{dln),^=^,,, (26) 

where the /th shell momentum fc; is given by 

I , ^ (/-!)(/ -2) m^" f-hl + h (m^r ,.7, 

^'-^° + — 16 — -W — ^2 kT ^ ^ 




Figure 2: The hypersurfaces k^ = ki on which the moments fi flow according to Eq. (p9D. 



and 
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Ki= ko+ 



{I - 1)(/ - 2)(/ - 3) m^" (/ - 1)(/ - 3)(/ - 5) (m^') 



+ 
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In particular k^ = k"^/ hq = k^^ is the semiclassical shell represented by the dispersion rela- 
tion ([T9|). One thus arrives at the following intuitive picture, illustrated in Fig. 2: taking 
moments of the QBE (|llD corresponds to discretizing momenta with a finite set of hypersur- 
faces as given by Eq. (p6D, the flow along which is described by the associated distribution 
function fi. 

After these considerations it is easy to show that the / th moment of QBE ( jTTD becomes 

2) 



u^dji + dji^, - [{d^kl)f,_, + ^^^^ 



mi)fi_, = Colli, 



(29) 



where we dropped the term k\\ ■ dg,.fi for simplicity (it can be always reinserted by the 
replacement udt — > ujdt + k\\ ■ 9^||). The first four equations (p9D , with / = 0, 1, 2, 3, together 
with the closure condition (p5|), form a closed set for functions /o, /i, /2 and f^. These 
equations are our main result. They can be used as a starting point for studying plasma 
dynamics in out-of-equilibrium situations when gradient approximation applies. 

Given the simple 4-point interaction term depicted in Fig. 1, the /th moment of the 



collision term appearing in ( [29| ) becomes 

Colh = -T>fi + r< {fi + k^,-'Ki) (30) 

where F^'and Tf are the spectral projections of the self-energies iS^ and iS^ defined by 

L Jp.k'^p' 

T< = ^ f {27:f5\k + p-k' -p')Up + l)U'fk', (31) 

^ Jp,k',p' 

where / = J duj(Pp\\/[(2'n-)^2psc] is the semiclassical equivalent of the Lorentz invariant three 
dimensional measure, and the a;-component of the 5-function should be taken as 6{6{kx)ki + 
^{Px)Psc—(^{K)Kc~^ip'x)Psc)- However, by the same approximation we made to arrive to Eqn. 
(pO|), we can set ki — > ksc in the delta function, which immediately leads to the identification 
rf'^ = r^'^ = r^'^. We win see below (c/. Eq. (^11) ) how, after a particular change of 
variables, F^'^ reduces to a more canonical form. In the more general case, when second 
order gradients are included. Colli depends on moments /;, /;_i, /;_2 and is a functional of /o, 
but not of higher order moments, so that including the gradient corrections to Colli would 
not spoil the closure property. 

4 Quantum coherence and semiclassical BE 

In order to separate the effects of quantum coherence in the dynamical equations (^91), it is 
convenient to define the following linear combinations of f]: 



f 


— ksc JO 




/qcl 


= fl-f 




/qc2 


= f 21 1^2 ~ /qcl " 


-/ 


/qc3 


= /s/^o ~ /qc2 " 


- /qcl 



/• (32) 

The distribution function / measures by definition population density of particles on the 
hypersurface kx = ko = ksc, while the densities f^d measure the correlations (quantum 
coherence) between neighboring hypersurfaces k^ = fcj-i and k^ = ki, where ki are defined 
in Eq. (|27|) . Note that in equilibrium in a uniform background all /jd vanish. 
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Multiplying Eq. ( pQ]) for /o by Kd^j one obtains 

dtf + —d.M + Aci) = -— / + — (/ + 1). (33) 

to to UJ 

This equation already resembles the standard Boltzmann equation, the main difference being 
coupling to an unknown function /qd . Dividing the /i-equation by oo and subtracting Eq. ( p3D 
one then obtains 

{dt + rqc)/qcl H -d^fqcl + — (/qcl + /qc2) H t^x/qc2 = Si, (34) 

to to to 

where the source 

si = dj + / 35) 

to to 

represents coherent mixing of / and /qd. One can obtain similar equations for the coherent 
quantum densities fqc2 and /qcs, but we will not present them explicitly here. The coherent 
density /qci is damped at the rate Fqc which reads 

r> _r< 

Tqc = . (36) 

to 



This is the out-of equilibrium generalization of the on-shell damping rate (c/. Eq. (|3lD and 
0). The coherence equations for /q^ are linear and hence can be quite easily solved, and 
the solution for f^ci inserted into Eq. (|33|). Requiring that none of /qd be sourced by the 
self-energy r< defines /qd uniquely. 

We now pause to discuss the validity of the gradient approximation. It is not hard to see 
that in the model under study the validity criteria reduce to the following conditions 

dtf « u;f, dj « kof, (37) 

which are the particular realization of ||0|| ^ 1 for the on-shell Boltzmann equation. 

In order to study how quantum coherence influences Eq. (0), we now make a simple 
estimate of /jd- Not far from equilibrium / can be approximated by its equilibrium form so 
that dxf ~ m'^/toLT, where L represents the scale on which m^ varies. To leading order in 
m^ the source si in Eq. (^) can be estimated as 

[LkoY ^ 
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Similar estimate holds for $2 and S3 in the equations for /qc2 and ^cs- If we want to include 
the higher order derivatives of k^c in Eq. ( p^ ) , this estimate implies that one cannot in general 
neglect /qd in Eq. (|3|). 

There is however a limit in which it is legitimate to neglect the part of /qd in Eq. ( P^D 
sourced by Si and still maintain higher order derivatives in ksc- Indeed, assuming efficient 
scattering, so that coherent quantum densities are strongly damped and the derivative terms 
in Eq. (^) can be neglected, we arrive at the following estimate of the coherent quantum 
densities sourced by Si. 



To check consistency of this estimate, note that in the efficient scattering limit one expects 
dxfqa ~ 1/-^) so that the spatial derivative term is suppressed by l/(rqc-L) in comparison 
to the Fqc^ci term. In the stationary frame the time derivative term is suppressed by an 
additional plasma velocity, which is often much smaller than unity. For example, the phase 
boundary speed at the electroweak phase transition is smaller than about 0.3 of the speed of 
light p. Upon inserting Eq. (^) into Eq. (|33|) , we see that the coherent quantum density /qd 
is suppressed by l/(rqcL) in comparison to the other third order gradient terms in Eq. (|33|) . 
We have thus shown that in the limit FqcL ^ 1 the coherent quantum density /qd in Eq. ( p3D 
can be neglected so that we arrive at the semiclassical Boltzmann equation 

dtf + ^dJ = Coll[fi (40) 

UJ 

where, at the leading order in gradients, the collision term is given in Eq. (^). Note that 
this equation still does not have the standard form since we have integrated over momentum 
kx and thus obtained a distribution function / = f{uj,k\^;t,x) = fcgc/o which is a function 
of energy u. To recover the semiclassical Boltzmann equation on the usual phase space we 
make the change of variables defined as follows 

ksc{uj;x) Cb{kx]x)' 

This choice is natural because it gives the local Lorentz covariant measure in the collision 
integral (^Tl). This change of variables in fact does not leave 5'^{k + p — k' — p') invariant. 

11 



The correction induced is, however, second order in derivatives of the distribution function 
moments in the coUision term, and hence it is beyond the approximation considered in this 
letter. Making use of dui -^ dCj{kx] x) = {dk^Lj)dk^, we then immediately obtain 

which specifies cD = tD(A;^; x). Consider now how the flow term transforms under this change 
of variables. We first have f{u] x, t) -^ f{kx; x, t), so that 



dx^ dx-{dk,uj)~^{dxuj)dk,, (43) 



and hence Eq. (^Of) becomes 

dtf + Vgdxf + FxdU = Coll[fl (44) 

with the following canonical velocity and force 

Vg = dk^u 

F, = -dxu. (45) 

This semiclassical Boltzmann equation is one of our main results. Note that "inversion" 
of the semiclassical dispersion relation k^^ = ksc{uj;x) as given in Eq. (^) is defined in a 
rather non-trivial way by invoking a natural transformation of the collision integral measure, 
resulting in a different form of the semiclassical energy hypersurface uj = u;(fc^; x) from what 
one might naively guess. This is the unique choice which, within the present approximation, 
renders the canonical form for the semiclassical velocity and force (^Sf). Recall again however, 
that in the present work the O^-term operating on the collision term and other analogous 
contributions were neglected. In a completely consistent computation these terms must 
be retained, and we are currently pursuing this computation in [Q. We finally emphasize 
that in case when the frequent scattering limit does not apply, one has to use more general 
equations (p9D with the closure relation (p5|). 
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5 Discussion 

We have generalized the semiclassical Boltzmann equation to include the effects of a slowly 
varying background to nontrivial order in gradients for a complex scalar field. Our motivation 
is baryon production at a first order electroweak phase transition, where this generalization 
is necessary to model CP violating effects which first appear in transport equations beyond 
leading order in gradient expansion. We found that consistent treatment of the system 
to nontrivial order in gradients requires the introduction of coherent quantum densities 
in addition to the usual distribution function, and derived the corresponding dynamical 
equations. We also showed that in the limit of frequent scatterings coherent quantum effects 
are suppressed, and the problem reduces to a single Boltzmann equation for the distribution 
function, containing a classical force which includes gradient corrections to nontrivial order. 
This is a particularly nice realization of the quantum-to-classical transition which occurs 
by the dissipative dynamics within the system itself. The dynamics is rendered dissipative 
by the weak coupling truncation of the Schwinger-Dyson equation at the order A^, and by 
the truncation of the gradient expansion at third order in gradients, resulting in localized 
dynamics in the Wigner representation and irreversibility. The rate at which the coherent 



quantum effects are dissipated is simply the out-of-equilibrium damping rate (pq). 

We now comment on how to relate our results to electroweak baryogenesis. The relevant 

CP violating effects emerge generally only at non-trivial order in gradients of the background, 

which is analogous to the case considered here. To make a connection with the computation 

of baryogenesis sources, it is instructive to integrate equation ([33| ) over uj and pjj to obtain a 

continuity equation for the density of particles which contains clearly separated semiclassical 

and quantum mechanical sources of the form 

. _ r du<fp\\ ksc . _ f du<fp\\ he . . 

•^^^ J (27r)3 cu^' ^""^ J (27r)3 cu ^'"='- ^ ^ 

This should be contrasted with the situation in literature 0, |[T0|, [|11|], where there is no 

agreement on whether the relevant CP violating sources come from classical or quantum 

currents. Moreover, in the frequent scattering limit, TqcL » 1, we find that the coherent 

quantum current jqc is suppressed, and the (dominant) semiclassical current then should be 

calculated making use of the semiclassical transport equation (0). This contains a classical 
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force analogous to the one first introduced as a source for baryogenesis in studies of two 



Higgs doublet models in Ref. 0, and then subsequently applied to the MSSM in Ref. [|T^ . 
Here we studied for simplicity a complex scalar field, whereas the physically relevant cases 
involve mixing scalar or fermion fields. However we believe that the novel features we found, 
being essentially a consequence of localization in space, are not specific to the complex scalar 
theory, but generic to all models. How this is realized in detail in other theories is under 
investigation |^, ||. 
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